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Abstract 

u' 

. In this work, we propose the design and analysis of a novel continuous robust controller for a class 

(-h of multi-input multi-output (MIMO) nonlinear uncertain systems. The systems under consideration 

contains unstructured uncertainties in their drift and input matrices. The proposed controller compensates 
the overall system uncertainties and achieves asymptotic tracking where only the sign of the leading 
minors of the input gain matrix is assumed to be known. A Lyapunov based argument backed up with 
an integral inequality is applied to prove stability of the proposed controller and asymptotic convergence 

m ■ 

of the error signals. Simulation results are presented in order to illustrate the viability of the proposed 



^1- 

o 



method. 



I. Introduction 

The problem of controller design for multi-input multi-output (MIMO) systems with uncer- 
tainties have been the research interest of many researchers. When the system under consideration 
is linear, solutions as early as flU are available. To name a few, in [1], assuming that the exact 
knowledge of high frequency input gain matrix is available, an adaptive controller was proposed. 
Ioannou and Sun [2 J used a less restrictive assumption on the gain matrix, the controller proposed 
required the existence of an auxiliary matrix which when pre-multiplied by the high frequency 
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gain matrix the resultant matrix would be positive definite and symmetric. In [3|, an adaptive 
controller for minimum phase systems with degree one has been proposed assuming that only 
the signs of the leading minors of the high frequency input gain matrix is known. For nonlinear 
MIMO uncertain systems the problem becomes more complicated. Therefore the solutions for 
only some special cases are considered in the literature. An adaptive backstepping method 
for strict feedback systems was utilized in [4] with the assumption that the input gain matrix 
premultiplying the control input is known. In [5], a general procedure for the design of switching 
adaptive controllers including feedback linearizable, and parametric-pure feedback systems has 
been proposed. An adaptive neural controller for MIMO systems with block-triangular form was 
proposed in [6]. Some more applications of controller design for MIMO uncertain systems like 
visual servoing, thermal management, aeroelasticity vibration suppression and ship control were 
presented in 0, M, M and [ED, respectively. 

Recently, |fTT|. |fT2|. [fT3l have proposed robust and adaptive type controllers for the MIMO 
nonlinear systems of the form 



x 



(n) - H (x, x, ■ ■ ■ , x {n -V) + G (x, x, ■ ■ ■ , x (n ~ 2) ) t (1) 



where (•)' denotes the i th time derivative, H (•) G R mxl and G (•) G R mxm are first order 
differentiable uncertain functions with G (■) being a real matrix with non-zero leading principal 
minors. Specifically in ifTD . authors have extended the work of llT4l by redesigning the controller 
of |[T5l removing an algebraic loop and potential singularity in their previous design and obtained 
a global Uniformly Ultimately Bounded (UUB) tracking error performance. In [fl2|. an adaptive 
controller that ensures asymptotic error tracking has been proposed. Recently, a continuous 
robust controller achieving semi-globally asymptotic tracking performance for uncertain MIMO 
systems of the form (H)) with two degrees of freedom was proposed in |fl3l . 

In this work, we consider a slightly broader class of uncertain MIMO nonlinear system than 
that of (0Q) which is of the form 

i (fl) = h(X) + g (X) t (2) 



1 T 



where X (t) = x T x T ■ ■ ■ (x (ri_1) ) r G R( mn ) xl is the combined state vector with 
x® (£) G R mxl i = 0, n, being states, h (■) G R mxl is an uncertain signal, g (■) G ]R mxm is a 
real matrix with non-zero leading principal minors, and r (t) G IR mxl is the control input. To 
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our best knowledge there are only few works on this model in the literature. Namely in [16], Xu 
and Ioannou considered the case where g (X) is either positive or negative definite, and designed 
a neural networks based adaptive controller that ensured local convergence of the tracking error 
to a residual set. While in IfTTl , Xian et al. considered the case where g (X) is positive definite 
and a robust controller containing the integral of the signum of the error term was designed to 
obtain semi-global asymptotic tracking. More recently in [TT8l . Chen et al. proposed a robust 
controller fused with a feedforward compensation term that ensured ultimate boundedness of the 
tracking error. Output feedback version of [|T8l were then proposed in |fl9l , and Il20l . 

In this paper, we have extended the full state feedback version of [18J to obtain asymptotic 
tracking as opposed to UUB. Apart from this, the proposed control design does not contain, 
though it is quite possible to design a similar one, an extra feedforward compensation term. An 
output feedback version can also be obtained similar to that of L19J, or [20J, however this result 
would possible lead an UUB error tracking performance. The main novelty of the method is 
the use of integral inequalities in conjunction with an initial Lyapunov based analysis to prove 
the boundedness of the error signals and then use the boundedness of the error terms in the 
final analysis to obtain the asymptotic tracking. The rest of the paper is organized as follows; 
Section II introduces the error system development while the full state controller development 
is introduced in Section III. Stability of the closed system under the proposed method and the 
numerical simulations are given in Sections IV and V, respectively. Finally, concluding remarks 
are presented in Section VI. 

II. Error System Development 

Based on the assumption that g (■) being a real matrix with non-zero leading principal minors, 
the following matrix decomposition is utilized [J3j , [1211 



where S (X) G M. mxm is a symmetric positive definite matrix, D e K mxm i s a diagonal matrix 
with entries being ±1, and U (X) G ]R mxm is a unity upper triangular matrix. Similar to [|T8l . we 
assume that D is available for control design. As we assumed that the leading principal minors 
of g (X) are non-zero, from (|2) it is straightforward to obtain 



g = S(X)DU(X) 



(3) 



T = g- 1 (x (n) - h) . 



(4) 
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Taking the time derivative of the system model in © and then substituting © results 

= h + gg- l (x^-h)+gf 

= ip + SDUr (5) 

where © was utilized and ip (X, x^) G R mxl is an auxiliary signal defined to have the following 
form 

( p = h + gg- 1 ( X W-h). 
Multiplying both sides of © with S* -1 (X) results in 

and after defining M (X) = S' 1 G R mxm and / (X, x^) = S^ip G M mxl , we obtain 

Mx (n+1) = j + DUf 

where M (X) is assumed to satisfy the following inequality 



(6) 



(7) 



(8) 



m Hxir < x t m (X) x <rh (X) ||xir v x e 



nmx 1 



(9) 



with m 6 R is a positive bounding constant, and fn (X) 6 K. is a positive, non-decreasing 
function. 

Our main control objective is to ensure that the system states would track a given smooth 
desired trajectory as closely as possible. In order to quantify the control objective, the output 
tracking error, ei (t) E R mxl , is defined as the difference between the reference and the actual 
system states as 

ei = x r — x (10) 
where x r (£) G R mxl is the reference trajectory satisfying following properties 

x r it) G C n , x® (t) G Coo , i — 0, 1, (n + 1) . (11) 



To ease the subsequent presentation, a combination of the reference trajectory and its time 



T • T 

3j 0u 



derivatives is also defined as X r (t) = 
design objective is to develop a robust control law that ensures 



(n-l) 



G R( mn ) xl . The control 



e« (t) 



— >■ as £ — >■ +oo, 



0, ...,n, while ensuring that all signals within the closed-loop system remain bounded. 
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In our controller development, we will also assume that the combined state vector X (t) is 
measurable. 

To facilitate the control design, auxiliary error signals, denoted by (i) G M mxl , i — 2, 3, n, 
are defined as follows 

e 2 = ei + ei (12) 
e 3 = e 2 + e 2 + ei (13) 



e n = e n -i + e n _! + e n _ 2 . (14) 

A general expression for (t), i = 2, 3, n in terms of e± (t) and its time derivatives can be 
obtained as 

i-1 

e i = ^ a i,i e l ] ( 15 ) 
j=0 

where ay G IR are known positive constants, generated via a Fibonacci number series. Our 
controller development also requires the definition of a filtered tracking error term, r (t) G M mxl , 
defined to have the following form 

r = e n + ae n (16) 



where a G R mxm is a constant positive definite, diagonal, gain matrix. After differentiating (1161 ) 
and premultiplying the resulting equation with M (•), the following expression can be derived 

Mr = M ^4 n+1) + an,A j+2) + j - / - DE/f (17) 

where d8]), (flOl) . (fT3T >. and the fact that a n) ( n _i) = 1 were utilized. Defining the auxiliary function, 

N (X,x( n \t) G M mxl , as follows 

N = M + a nA +2) + «e„ j ~f + e n + X -Mr (18) 

the expression in (fTTb can be reformulated to have the following form 

Mr = —Mr - e n - DUf + N. (19) 



Furthermore, the filtered tracking error dynamics in (1191 can be rearranged as 

Mf = --Mr-e n -D(U-I m )i-Di + N + N (20) 

January 24, 2013 DRAFT 



6 



where we added and subtracted Df (t) to the right-hand side, I m G R mxm is the standard 
identity matrix, and N (t), N (i) G W mxl are auxiliary signals defined as follows 

N = N-N. (22) 

The main idea behind adding and subtracting Df (£) term to the right-hand side of (|20l ). is to 
make use of the fact that U (•) is unity upper triangular, and thus (U — I m ) is strictly upper 
triangular. This property will later be utilized in the stability analysis. 

III. Controller Formulation 



Based on the open-loop error system in (1201 ) and the subsequent stability analysis, the control 
input, r (£), is designed to have the following form 



t 

e„ (t) - e n (t ) + a I e n (a) da 



to 



DU (23) 



t = DK 

where the auxiliary signal IT (t) e R mxl is generated according to 

n = CSgn(e n ) j n(t ) =0 mxl . (24) 



In (1231) and (I24|) . K, C G R mxm are constant, diagonal, positive definite, gain matrices, mx i G 
R mxl is a vector of zeros and Sgn(-) G R mxl is the vector signum function. Based on the 
structures of (|23l and (|24l . the following expression is obtained for the time derivative of the 
control input 

t = DKr + DCSgn (e n ) (25) 

where (fT6l ) was utilized. The control gain is chosen as K = I m + k p I m +dia.g {kd,i, ^.(m-i), 0} 
where k p , k^% G R are constant, positive, control gains. Finally, after substituting (1251) into (|20l) . 
the following closed-loop error system for r (i) is obtained 

Mr = -^Mr - e n - iTr + N + N - D (U - I m ) DKr - DUDCSgn (e n ) (26) 

where the fact that DD = I m was utilized. 

Before proceeding with the stability analysis, we would like to draw attention to the last two 
terms of (|26l ) which we will investigate separately in the next two subsections: 
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A. The "D(U - I m )DKr" Term 

Note that, after utilizing the fact that (U — I m ) being strictly upper triangular, we can rewrite 
the term D (U — I m ) DKr as follows 

A + $ 



D(U-I m )DKr 







where A(t), $ (t) G IR (m 1)xl are auxiliary signals with their entries A, (t), $j (t) G 
1, (m — 1), being defined as 

Aj = di 2J djkjUijTj 
j=i+i 
m 



j=i+i 



with C/jj (X r ), [/jj (t) G R are defined as 



U ■ ■ = U- ■ 



By 



i,j|X=X r 



(27) 



(28) 
(29) 



(30) 
(31) 



where £/y (X) G 1 are the entries of U (X). Notice from (ITTT) that the last entry of the term 
D {U — I m ) DKr is equal to 0, and its z-th entry depends on the [i + l)-th to m-th entries of 
the control gain matrix K. 



B. The "DUDCSgn (e n ) " Term 

We can rewrite the DUDCSgn (e n ) term as 







+ e 



DUDCSgn (e n ) = 

where ^ (t) G ]R( m_1 ) xl and (t) G R mxl are auxiliary signals defined as 

= D(U-U)DCSgn(e n ) 
6 = DUDCSgn (e n ) 







(32) 



(33) 
(34) 
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where U (X r ) = U\x=x r £ M mxm is a function of reference trajectory and its time derivatives, 
and (t) £ K, i = 1, (m — 1) and 0j (t) e R, i = 1, m, are defined as 



= di ^ d jCjUijSgn(e n j] 



(35) 



©i = ^ djCjUjjsgn (e n j 



(36) 

Remark 1: The Mean Value Theorem [|22| can be utilized to develop the following upper 
bounds 



N(-] 



< 



PjvlFlU W z \ 



< /MINI) IN 



(37) 
(38) 



where pjy (•)> Pij (•) eKare non-negative, globally invertible, non-decreasing functions of their 
arguments, and z (t) E R^ n+1 ^ xl is defined by 



T T 
e l e 2 



T T 



(39) 



It can be seen from CCD), (fl8j (|2TT) that N (t) and L/^- (t) are bounded in the sense that 0/7) 

\Ni(t)\ < Cn, (40) 

\UiAt)\ < (41) 
where (j^., e R are positive bounding constants. Based on (1281) . (|29l , (|35l) . (|36l) . following 
upper bounds can be obtained 

rn 

\K\ < Yl k JpiA\\ z \\)\\4\rj\<PAA\\4)\\4 



j=i+i 

m 



m 



j=i+l 
m 



I ©,| < E^^Ce, 



(42) 
(43) 
(44) 
(45) 



where (|37T) -(|4TT) were utilized. From (I43T ). it is easy to see that ||6|| < Ce is satisfied for some 
positive bounding constant (e 6 R, and from (142]) — (144]) . we have 

1^1 + 1^1 + 1^1 < a (ikii)ikii (46) 
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where (||z||) G R z = 0, 1, (m — 1), are non-negative, globally invertible, non-decreasing 
functions satisfying 

Pa,. + PVi + C*» < Pi- (47) 

Remark 2: Notice that, as a result of the fact that U (t) being unity upper triangular, 6 (t) in 
( 1341 ) can be rewritten as 

e=(/ w + fi)CSgn(e n ) (48) 

where Q(t) = D (U — I m ) D G R mxm is a strictly upper triangular matrix. Since it is a function 
of the reference trajectory and its time derivatives, its entries, denoted by fiy (£) G R, are 
bounded in the sense that 

|a,il<C^ (49) 
where £ n G R are positive bounding constants. 

At this point, we are now ready to continue with the stability analysis of the proposed robust 
controller. 



IV. Stability Analysis 

This section, via an initial Lyapunov based analysis we will first prove the boundedness of the 
error signals under the closed-loop operation. Using this result we will then present a lemma 
and obtain an upper bound for the integral of the absolute values of the entries of e n . This upper 
bound will later be utilized in another lemma to prove the non-negativity of a Lyapunov-like 
function that will be used in our final analysis which proves asymptotic stability of the overall 
closed-loop system. 

Theorem 1: (Boundedness proof) For the uncertain MEMO system of 0, the controller in 
(1231) and (T24l) guarantee the boundedness of the error signals (flOl) . (fl2l)-(fT4"l) and (fT6|) provided 
that the control gains kd,i and k p are chosen large enough compared to the initial conditions of 
the system and the following condition is satisfied 

A min (a) > - (50) 

where the notation A min (a) denotes the minimum eigenvalue of the gain matrix a, previously 
defined in (fl6i 

Proof: See Appendix |Al ■ 
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Lemma 1: Provided that e n (t) and e n (t) are bounded, the following expression for the upper 
bound of the integral of the absolute value of the i-th entry of e n (t) i = 1, • • • , m can be 
obtained 

t t 
J |e„,i(a)|rf(j < 7i + 72 J \e Uji (o)\ da + \e n>i \ (51) 

to to 

where 71, j 2 G M are some positive bounding constants. 

Proof: The proof is very similar to that of the one given in [23], however for the complete- 
ness of the presentation we have included it in Appendix El ■ 
Lemma 2: Consider the term 

L = r T (iV — (J m + ft) CSgn (e n )) (52) 

where ft (t) defined in (|48l ) is strictly upper triangular and also is a function of desired trajectory 
and its time derivatives; thus it is bounded. Provided that the entries of the control gain C are 
chosen to satisfy 

C m > <N m (l + ^-) (53) 

C l > |u + J2 Cni^iJ (l + . * = 1> (™ " !) (54) 
then it can be concluded that 

jL(a)da<( L (55) 

to 

where (j, G M is a positive bounding constant defined as 

m— 1 m mm 

(l = 7i^E Cn.jC, + 71 E + E ^ (*>) I • (56) 

i=l j=i+l i=l i=l 

Proof: See Appendix O ■ 
Theorem 2: (Asymptotic convergence proof) Given the uncertain MIMO nonlinear system 
of the form (0), the continuous robust controller of (|23l and ([24]) ensures that all closed-loop 
signals remain bounded and the tracking error signals converges to zero asymptotically in the 
sense that 

->• as t ->■ +00 , Vi = 0, • • • , n 
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provided that a is chosen to satisfy (|50l ), the entries of C are chosen to satisfy (|53l and (|541 ), 
and A; rf j and fc p are chosen large enough. 

Proof: See Appendix iDl ■ 

V. Numerical Results 

In order to substantiate the theoretical results, numerical analysis has been carried out. Similar 
to ffi~8l the performance and liability of the proposed nonlinear robust controller has been tested 
on a two-link robot manipulator system with coupling between the two links taken from [|24|. 
The equations of motion are given as lfi~8l 

-hq 2 -h(qi + q 2 ) 
-h qi 

where qi (t), q 2 (t) G R denote the positions of the joint angles, and H n , H u , H 22 and h are 
explicitly defined as 









H\ 2 




9i 












+ 






H\ 2 


H 22 




fa 







9i 




92 



(57) 



H\ 2 
H 22 
h 



ai + 2a 3 cos q 2 + 2a^ sin q 2 
a 2 + a 3 cos q 2 + a 4 sin q 2 

a3sing 2 — a4sing 2 - 

where a\ = 4.42, a 2 = 0.97, a 3 = 1.04 and a 4 = 0.6. The term 
as 

T* 

/3(9i,92) 



TV To 



To 



1 1 

1 



n 

T 2 



(58) 
(59) 
(60) 
(61) 

in (1571) is obtained 
(62) 



where ri (t) and r 2 (t) are the control inputs, and j3 = HnH 22 — Hf 2 e R. The robot manipulator's 
initial positions have been set to q (0) = [ 10 10 ] T (deg) and g(0) = [ ] T (deg.s _1 ). The 
control objective is to make q± (t) and q 2 (t) follow a sinusoidal desired trajectory chosen as 

q d {t) = (1 -exp(-0.3t 3 ))[ 30sin(t) 45sin(t) ] r (deg). 



After a rough tuning process the control gains are selected as follows 

a = diag{l, 5}, K = diag{l75, 125} and C = 5. 
where diag {•} is used to represent a diagonal matrix. 
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Remark 3: Though in the simulations the tuning process did not require the separation of 
gains k dji and k p from K on an actual system this might be needed. In these cases, adjusting 
control gains might become a bit tricky. We advice the users to adjust the entries of C first, then 
based on C adjust k p , and finally work on k d based on k p . 

The link position tracking error is depicted in Figure [TJ while the control input is shown in 
Figure |2l Simulation results confirm that the proposed controller meets the position tracking 
objective and asymptotic tracking is achieved. 
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Fig. 1. Link Tracking Errors 





VI. Conclusion 

In this work, for a class of uncertain MIMO systems having non-zero leading principle 
minors in their input gain matrix, a continuous nonlinear robust controller have been proposed. 
Via the use of Lyapunov based arguments in conjunction with an integral inequality we were 
able to obtain semi-global asymptotic tracking. The proposed controller can compensate the 
uncertainties though the system under a smoothness assumption on the system uncertainties. 
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Fig. 2. Control Inputs 

Simulation results are also included to illustrate the viability and the performance of the proposed 
method. 
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Appendix A 
Proof of Theorem d] 

Proof: The non-negative function V\ (z) E R is defined as 

1 n 



efet + -r T Mr. 



(65) 



i=l 



(66) 



By utilizing ©, (1651) can be bounded in the following manner 

1 2 1 _ 2 

— min {1, m} ||z|| < Vi (2) < — max {1, m (||z||)} ||^|| , 

where z(t) was defined in (1391) , and the terms m , m(\\z\\) were defined in (|9]). Taking the time 
derivative of (1651 yields 

1 

Vi = efe* + r T Mr + -r T Mr. 



(67) 



8=1 



The first term in the above expression can be written as follows 

n 

2J ef e, = (e 2 - ei) + (e 3 - e 2 - 



i=i 



+ej (e 4 - e 3 - e 2 ) + ... 



+e^_! (e n - e n _i - e n _ 2 ) + (r - ae n ) 



n-l 



ET T T 1 

e i + e n _ 1 e n + e n r — e r 



(68) 



i=i 



where (IT2l)-(IT4l). (1161) were utilized. Substituting (I26l)-(l29l). (1321)— (1341) and ([68]) into d67j) results 
in 

n-l 

2=1 

,T 



: /•' ( -^Mr - e n - Kr + N + N 



A + $ 


which, after substituting the control gain matrix, can be rewritten as 



T 




— r 










r T 6 + -r T Mr 
2 



(69) 



n-l 



Ti 



T 



i=l 



ET T T 

6j 6j -|- 6 n _-^6 n 6 n OC n 

m— 1 



+ r N — k p r T r 
+r T N - r T Q. 



+ 



m— 1 



i=l 



i=l 



(70) 
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After completing the squares in bracketed terms, bounding N (£) and (t) with constants, and 
utilizing e^_ x e n < 1/2 ||e n _i|| 2 + 1/2 ||e n || 2 , we obtain 

V\ < ~ ^ |N| 2 ~ \ IK-if ~ y A min (a) - ||e n || 2 - r T r 



i=l 

n 2 /II _ 1 1 \ to— 1 



, ffrira „ „2 , a- I, ,.2,11 I, > ,n I, t ~ n 



which can then be rearranged as 



^<_f Al _»_£4P)|| 2f + ^ (72) 

where A x = min {|, A min (a) — |, 1 — 5 G R is a positive bounding constant, e = Cat + Ce> 
and ||r|| e < i IM| 2 + fc 2 was utilized. When the controller gains kd,i and k p are selected large 
enough (compared to the initial conditions of z (t)), and utilizing (|66l) , the following inequality 
can be obtained 

Vi < -p x V x + 6e 2 (73) 

where j3 x E R is a positive constant. From (l65"b . and (|73T) . we can conclude that V x (t), therefore 
ei (t) for i — 1, n and r (t) are uniformly ultimately bounded. ■ 

Appendix B 
Proof of Lemma CD 

Proof: First, we note that if e n> j(t) = on some interval, then e n i (t) = on the same 
interval, and the inequality (|5"TT ) yields this qualification. Therefore, without loss of generality, 
we assume that e n i (i) is absolutely greater than zero on the interval of [t ,t\. Let T e [to,t) 
be the last instant of time when e nji (t) changes sign. Then, on the interval [T,t], e nii (t) has a 
constant sign, hence 



T 



e nA (a)\da 



e„ ,■ (a) da 



T 



\e n>i (t)-e n;i (T)\. (74) 



From the boundedness of the function e n i (t), it follows that there exist a constant 7 > such 

that \e ni (t)\ < 7, therefore 

rT 

\e n>i {a)\da<j{T-t ). (75) 



to 
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On the other hand, we obtain the following equality from the application of the Mean Value 
Theorem [ |22l 

T 

\e n>i (a) | da = e n>u (T - t ) (76) 

to 

where e n ^ is some intermediate value of |e n i (7j)| on the interval [£o>^1- By assumption, e n it is 
bounded away from zero. Therefore we can conclude as follows by using inequality (f75]) and 
equality (|76|) 

/ \e n ,i (cr)| da < 72 / | e n>i (a) \ da (77) 

J to J to 

where 72 = 7/e^. Combining the relationships in (|74T ) and (1771 , we can write 

/ \e n ,i (a)\da < \e n:i (t)\ + \e nA (T)| + 72 / |e„,» (cr) | da (78) 

■/ to J to 

which yields the inequality (|5TT) with definition 7 X = sup |e n i (T)|. ■ 

Appendix C 
Proof of Lemma [2] 



Proof: We start our analysis by integrating (|52t in time from t to t 

t t 
J L(a)da = J el (a) a T (N (a) - CSgn (e n (a))) da 



to to 

t 



ei (a) a 1 ft (a) CSgn (e n (a)) rfa 



to 



t t 
+ 1 e T n (a)N(a) da - J el (a) ft (a) CSgn (e n (a)) da 



to t a 
t 



J el (a) CSgn (e n (a)) da (79) 



to 



where (fT6l) was utilized. To ease the presentation, we will consider each term on the right-hand 
side of (|79l) separately: 
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The First Term: 



„ » m 

j e n (?) ® T (N ( a ) - CSgn ( e n (<t))) da = Yl aie ' n ^ ^ ~ ^S 11 ( en >* rf(J 



to to 



* =1 



(80) 



The Second Term: 



t 1 

, 171 — 1 



II L ± lit 

e T n (a) a T tt (a) CSgn (e n (cr)) da = - / c^e^ (<r) ^ (cr) sgn (e nj - (ex)) dcr 

1 t 

m— 1 m „ 

< EE ^ C ^h,j \e n A°)\d° (81) 



i=i j=i+i t0 



77ze TTzzrd Term: 



t t 

m 



/p III: 
e T n (a) N (a) da = (a) N f (a) da 

r i=i 



to to 



m „ 

- E^ / |en,i(^)|^ 
1=1 / 

^ E Ctv, ( 7i + 72 y |e„,i (a) | da + |e„, 



(82) 



to 
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The Fourth Term: 



* m-l 



/n 1 1 V A (ft. 
(a) (a) CSgn (e n (a)) dtr = ~ / E ^ ( a ) E ^'^J (°") s § n ( e "J (°")) dcr 



to to 



i=l j=i+l 



1 * 
m—l m p 

< EE J (*) 



da (83) 



•=i j=i+i 

m— 1 m 



i=l j=i+l 

The Fifth Term: 



< E E C iC^ij |7i + 72 J \e n ,i (a)\da + \e n ^ 



to 



t 



to to 



/p lit! 
el (a) CSgn (e n (<r)) tZtr = - / E (°0 s S n ( e rM (°0) rf(T 

to 1=1 

t 

= "E^ / sgn(e nji (<r))d(e n)i ) 
i=i / 

to 

= -£cw 

i=1 to 

m m 

= -E^K^)i+E c *i e ^(*°)i- < 84) 



i=\ i=l 
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It is noted that, the result of Lemma \T\ was utilized to obtain (1821 and (f83T ) . After combining the 
upper bounds in (l8~0l)-(l8~4T). we obtain 



/ L{a)da < + E <*u C i) - C * I KiWI^ 



72 



a, 



fa - ) I da 



to 



+ (Cat™ — C m ) \e n<r . 



m—l 



(85) 



i=i \ j=t+i / 

m—l m mm 

+^ 2 S + 7i£)u + E Ci i e ™<* ( to )i • 

8=1 j'=i+l i=l i=l 

Based on (1851) . we first choose C m to satisfy (1531) to make second and third expressions on the 
right-hand side negative, we next choose Cj starting from (m — l) with a decreasing order to 
satisfy (l54l) to make first and fourth expressions on the right-hand side negative, and finally, we 
utilized the definition of ( L in (l56l) to obtain (1551) . thus completing the proof of Lemma HI ■ 



Appendix D 
Proof of Theorem [2] 

Proof: Let the auxiliary function P (i) e R be defined as follows 

P = Cl- fL(a)da. 

Jt 



(86) 



where the terms £l and L(i) were defined in (|56l) and (|52l) . respectively, when the entries of the 
control gain matrix C are chosen to satisfy (1531) and (|54l) . from the proof of Lemma |2] given in 
Appendix O we can conclude that P (t) is non-negative. At this stage, consider the Lyapunov 
function, denoted by V (s, t) e K, defined as follows 



V = V x + P 



where s (t) G l [( " +1)m+1]xl is defined as 



A 

s = 



1 T 



P 



(87) 



(88) 
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By utilizing ©, (f8VT ) can be upper and lower bounded in the following form 

W 1 {s) <V(s,t) < W 2 (s) 
where Wi(s), W 2 (s) G M. are defined as 

W l ±\ 2 \\s\\\W 2 ±\ 3 (\\s\\)\\s\\ 2 



(89) 



(90) 



with A2 = | min {1, m} and A3 = max {l, \m (||z||) } . 

Taking the time derivative of V (t), utilizing the time derivative of (|55l ). canceling common 
terms and following similar steps to that of proof of Theorem Q] yields 



n-1 



V 



ET T T T 

&i &i ^n—l^n ^ n Q!6 n T T 



i=l 



+ 



r T N — k p r T r 



m—l 



TO— 1 



8=1 



i=l 



(91) 



which can be rearranged to have the following form 



n-2 



E 

i=l 



1 



V < - > lk:l| 2 - - ||e n _i|| 2 



Amin (a) - - ) ||e n || 2 - r T r 



m—l 



MI 2 + £ 



< 



A, 



i=l 

m-1 9 



E 



pKWA 



4k p Ak d)i 



\\4 



(92) 



(93) 



where A 4 = min {-, A min (a) — |}. When the controller gains k p and k d ,i for i = 1, 2, (n — 1) 
are selected large enough such that the regions defined by V z = {z : ||2|| < 7?.} and V s = 
{s : < 72.} with 72 defined as 



1-/3 



-1 



m 



2\ h 



1-/3 



for i = 1,2,..., (m-1) (94) 



are non-empty, from (1931 and the definition of s, one can restate 



V < -0\\z\ 



-W(s),Vs G X? s 



(95) 



where /3 G K is a positive constant that satisfies < /3 < 1. From the definition of (|87l ) and 
(1931 ), it is obvious that and V (t) G Coc, also from the proof of Theorem [1] and outcome of 
standart linear analysis methods, we can conclude that all signal in the closed-loop error system 
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are bounded and furthermore, from the boundedness of W(s), we can state W (s) is uniformly 
continuous. 

Based on the definition of V s , another region, S, can be defined in the following form 



S= |^P s :^ 2 ( S )<A 3 (p~ 1 (2 V / ^)y| 
n Is eV s :W 2 (s) < A 3 (p? (2yjk^i=f)^ 

ni 



(96) 



1-/3 




A direct application of Theorem 8.4 in [22] can be used to prove that \\z (t)\\ — > as t — > +oo 
Vs (to) £ <5. Based on the definition of z(t), it is easy to show that ||ej(t)|| , \\r (t)|| — > as 
t — )■ +oo Vs (t ) £ 5, i = 1,2, ...,n. From (fT6l . it is clear that ||e n (t)|| — > as t — > +oo 
Vs (t ) £ <S- By utilizing (fl"5T ) recursively, it can be proven that (t) — )■ as t — > +oo, 
i — 1,2, ...,n Vs (t ) £ 5. Note that the region of attraction can be made arbitrarily large to 
include any initial conditions by choosing the controller gains k p and k d:i for i = 1,2, (m — 1). 
This fact implies that the stability result obtained by proposed method is semi-global. ■ 
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